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Abstract 

Background: Advanced paternal age is robustly associated with several human neuropsychiatric disorders, 
particularly autism. The precise mechanism(s) mediating the paternal age effect are not known, but they are 
thought to involve the accumulation of de novo (epi)genomic alterations. In this study we investigate differences 
in the frontal cortex transcriptome in a mouse model of advanced paternal age. 

Findings: Transcriptomic profiling was undertaken for medial prefrontal cortex tissue dissected from the male 
offspring of young fathers (2 month old, 4 sires, n= 16 offspring) and old fathers (10 month old, 6 sires, n= 16 
offspring) in a mouse model of advancing paternal age. We found a number of differentially expressed genes in 
the offspring of older fathers, many previously implicated in the aetiology of autism. Pathway analysis highlighted 
significant enrichment for changes in functional networks involved in inflammation and inflammatory disease, 
which are also implicated in autism. 

Conclusions: We observed widespread alterations to the transcriptome associated with advanced paternal age 
with an enrichment of genes associated with inflammation, an interesting observation given previous evidence 
linking the immune system to several neuropsychiatric disorders including autism. 
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Findings 

Background 

Advanced paternal age is robustly associated with several 
neuropsychiatric disorders, most notably autism [1]. The 
mechanism(s) underlying this paternal age effect are not 
known, although the accumulation of de novo (epi) gen- 
etic changes over time in the male germ line are thought 
to be important. Evidence to support this hypothesis 
comes from preliminary data showing a higher rate of de 
novo mutations [2] and altered DNA methylation at 
birth [3] with increasing paternal age in humans, and in- 
creased de novo copy number variation (CNV) [4] and 
altered brain DNA methylation at imprinted loci in ro- 
dent models of advanced paternal age [5]. 

Gene expression changes associated with aging in 
the brain have been investigated largely in terms of 
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neurodegeneration and dementia [6], although aging per 
se is associated with widespread transcriptional changes 
[7]. In a study of aging using mice, for example, brain 
samples from the cerebellum and neocortex had in- 
creased expression of genes related to inflammatory and 
stress responses and decreased expression of genes asso- 
ciated with growth and trophic factors, protein turnover, 
DNA synthesis and repair, and neurotransmission [8]. 
Age-associated changes have been observed in numerous 
tissue and cell types, including the germ line. Of relevance 
to paternal age, an expression study of rats of different 
ages identified over 2,800 loci that are differentially 
expressed in spermatocytes from older males (18 months) 
compared to spermatocytes from young males (4 months); 
of note, many genes associated with base excision repair, 
nucleotide excision repair, mismatch repair and double 
strand break repair were altered in spermatocytes from 
older males [9]. In an analysis of RNA from the spermato- 
gonial stem cells of mice of four different ages (6 days, 
21 days, 60 days and 8 months) using microarrays, 2,819 
genes had differential expression between the age groups 
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{P < 0.05, fold change >2) including genes previously iden- 
tified in gene expression studies of aging in stem cells [10]. 
Pathway analysis of these genes highlighted an enrichment 
of genes involved in DNA repair and oxidative stress, 
which is interesting given the known increase in DNA 
damage in the spermatozoa of older males. To date, only 
one study has looked at gene expression changes associ- 
ated with paternal age in humans [11]. In this study of pa- 
ternal age and autism, a decrease in the overall variance of 
gene expression was observed in the offspring of older fa- 
thers, in addition to a downregulation of genes involved in 
gene transcription [11]. To date, no research has examined 
transcriptomic changes in the brain in the context of ad- 
vanced paternal age. 

Methods 

C57BL/6J mice were bred and maintained in the Bio- 
logical Services Unit at the Institute of Psychiatry, Kings 



College London, using stocks purchased from Charles 
River Laboratories. All animal experiments received the 
approval of the local ethical review panel of King's 
College London and were performed in compliance with 
the UK's Animals (Scientific Procedures) Act 1986. The 
work was carried out under licence (PPL 70/7184) and 
all efforts were made to minimize animal suffering and 
to reduce the number of animals used. We used tissue 
collected from an experimentally controlled rodent 
model of advancing paternal age (described in detail pre- 
viously [5]) to examine whether the offspring of older fa- 
thers have altered levels of gene expression. Male 
offspring of 'young' fathers (2 month old, n = 16, four 
sires with one litter each) and 'old' fathers (10 month 
old, n = 16, six sires with one litter each), with at least 
two individuals selected from each family, were used in 
this study (see Additional file 1 for an overview of the 
samples used). All dissections were performed by a 
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Figure 1 Top-ranked differentially expressed transcripts between the offspring of young fathers and offspring of old fathers. Shown for 
each probe (associated with the genes AA408296, Mud 5 and Lta) are (A) group-level differences and (B) data split by family. The y-axis shows the 
log2 transformed signal intensity. 
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single individual to ensure topographical similarity be- 
tween samples. Briefly, brains were removed from the 
skull and placed dorsal side down on a wetted filter 
paper on a petri dish kept on ice. The cerebral halves 
were opened out from the midline, after cutting through 
the corpus callosum. Approximately 3 mm^ of tissue 
was cut from the anterior part of the frontal lobes (from 
bregma 2.46 mm to 1.34 mm [12]), mainly containing 
the medial prefrontal cortex including some prelimbic 
cortex, infralimbic cortex, cingulate cortex and motor 
cortex [13]. 



High quality RNA (average RNA integrity number 
(RIN) 8.4) was isolated using the Qiagen AllPrep DNA/ 
RNA Micro Kit (Qiagen, Manchester, UK). Gene expres- 
sion was quantified using the Illumina Mouse Ref8 V2 
array (Illumina, San Diego, CA, USA), which targets ap- 
proximately 25,600 transcripts and over 19,100 unique 
genes. After stringent quality control and preprocessing, 
raw expression data were batch-corrected using the Com- 
Bat package [14] before being normalized and analysed 
using lumi [15] within the R statistical environment. All 
analysis scripts are available on request from the authors. 



Table 1 Probes showing gene expression differences associated with paternal age (P< 0.001) 



Transcript 


Probe ID 


Expression difference^ 


P value 


Mean (SEM) young 


Mean (SEM) old 


AA408296 


4150356 


T 


3.23 X 10"*^ 


6.62 (0.01) 


6.71 (0.01) 


MuclS 


2120176 


I 


7.45 X 10"*^ 


6.5 (0.01) 


6.42 (0.01) 


Lta 


2900474 


T 


7.46 X 1 0"*^ 


6.52 (0.01) 


5.61 (0.01) 


Tsip 


6960520 


T 


3.02 X 10"^ 


641 (0.01) 


5.48 (0.01) 


Zfp35 


580168 


I 


1.30X 10"^ 


6.99 (0.02) 


5.89 (0.01) 


Gstml 


4540458 


T 


1.31 X 10"^ 


6.57 (0.02) 


5.66 (0.01) 


Rdh20 


6840020 


T 


1.38X 10"* 


6.48 (0.01) 


5.56 (0.01) 


Ccdc90b 


6100215 


T 


1.91 X 10""* 


7.1 (0.02) 


7.22 (0.02) 


Sdhc 


3290630 


i 


2.10 X lO"** 


6.48 (0.01) 


6.4 (0.02) 


AtpSsI 


5270131 


i 


2.97 X 1 0"** 


7.76 (0.02) 


7.61 (0.03) 


Oifrl87 


380370 


i 


3.07 X lO"** 


6.59 (0.01) 


5.51 (0.01) 


Cd44 


4880138 


I 


3.29 X 10"* 


6.64 (0.01) 


5.57 (0.01) 


Fryl 


510059 


T 


3.58 X 10* 


6.57 (0.01) 


5.64 (0.01) 


Accnl 


7330561 


T 


3.59 X 10""* 


6.69 (0.01) 


5.76 (0.01) 


Cdl77 


670452 


i 


3.71 X 10"* 


6.51 (0.01) 


5.43 (0.01) 


Tspan3 


3130497 


T 


3.96 X 1 0"* 


10.33 (0.04) 


10.49 (0.03) 


Olfr92 


2710201 


T 


4.29 X 1 0"* 


6.47 (0.01) 


5.53 (0.01) 


Shh 


2260450 


T 


4.37 X 10"* 


6.71 (0.02) 


5.79 (0.01) 


E130303B06RIK 


6250669 


I 


4.38 X 1 0"* 


6.67 (0.01) 


5.59 (0.01) 


ArhgapS 


7550333 


I 


4.43 X 10"* 


6.45 (0.01) 


5.39 (0.01) 


Zc3hl8 


7380382 


T 


5.21 X 10"* 


7 (0.01) 


7.08 (0.02) 


AtgS 


3850450 


i 


6.10 X 10"* 


7.98 (0.02) 


7.9 (0.01) 


Nsfllc 


7610433 


I 


7.45 X 10"* 


7.69 (0.02) 


7.55 (0.03) 


Earl2 


4540543 


I 


7.47 X 10"* 


6.52 (0.01) 


5.46 (0.01) 


Olfrl366 


6860093 


I 


7.69 X 10"* 


6.49 (0.02) 


541 (0.01) 


Tnfsfl 1 


2480255 


T 


8.20 X 10"* 


6.5 (0.01) 


5.57 (0.01) 


Lcorl 


2470343 


i 


8.27 X 1 0"* 


6.56 (0.01) 


5.49 (0.01) 


Dusll 


2370138 


T 


9.30 X 1 0"* 


7.13 (0.01) 


7.21 (0.02) 


Larp4 


6400025 


i 


9.59 X 1 0"* 


6.47 (0.01) 


541 (0.01) 


Algl2 


3850500 


T 


9.79 X 10"* 


6.56 (0.01) 


5.65 (0.02) 


lerSI 


4900053 


I 


9.82 X 1 0"* 


7.6 (0.02) 


747 (0.03) 


Ccl26 


4560164 


I 


9.94 X 1 0"* 


6.49 (0.01) 


541 (0.02) 


^Arrows denote direction of change in 


the offspring of older father compared to the 


offspring of younger fath 


ers. 





SEM, standard error of the mean. 
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Results and discussion 

We found numerous differences in gene expression as- 
sociated with advanced paternal age. Most notably, 
probes associated with three transcripts were found to 
be differentially expressed in the offspring of older fa- 
thers compared to the offspring of younger fathers with 
a false discovery rate (FDR) < 0.1: AA408296 (also known 
as Diexf) {P = 3.23 x 10'^, FDR = 0.06) and Lta {P = 
3.23 X 10", FDR = 0.06) had elevated expression, whilst 
MuclS {P = 7.46 X 10 ^ FDR = 0.06) had reduced expres- 
sion (see Figure lA). Although the absolute differences 
in gene expression are small and future work is needed 
to ascertain their functional significance, for each of 
these differentially expressed transcripts, consistent pat- 
terns of altered expression were observed across each of 
the old-father families (Figure IB and Additional file 2). 
The top-ranked differentially expressed transcripts (P < 
0.001) are listed in Table 1. Strikingly, almost a quarter 
of these differentially expressed loci have been previously 
implicated in the aetiology of autism, including Gstml 
(deletion in GSTMl observed in autism case-parent 
trios) [16], Ccdc90b (missense mutation associated with 
autism in a study of cases and their parents) [17], Cd44 
(differential expression in autism discordant monozy- 
gotic twins) [18], Accnl (multiple SNPs associated with 
autism in a familial study) [19], Shh (increased serum 
SHH expression in autistic patients compared to non- 
autistic controls) [20], Dusll {de novo missense mutation 
observed in autism families) [21] and lerSl (increased ex- 
pression in lymphoblast cells from autism patients in 
discordant sibling pairs) [22]. 

Ingenuity pathway analysis (IPA) [23] and the Database 
for Annotation, Visualization and Integrated Discovery 
(DAVID) [24,25] were used to identify gene pathways and 
functions enriched amongst transcripts differentially 
expressed in the offspring of old fathers. Of note, IPA iden- 
tified a significant enrichment for functional pathways in- 
volved in inflammation and inflammatory disease. Table 2 
shows the top-ranked biological functions enriched in 
probes differentially expressed (P < 0.01) in the frontal cor- 
tex in the offspring of old fathers. The top-ranked biological 
network regulates immune cell trafficking and cell-to-cell 
signalling (Figure 2). Furthermore, the most significantiy as- 
sociated gene ontology term identified by DAVID was im- 
mune response (14 genes, P = 2.87 x lO '*) and the second 
most associated cluster (after cytokine activity) included 
defence response, response to wounding, inflammatory re- 
sponse and acute inflammatory response. This is interesting 
given mounting evidence linking the immune system to 
several neuropsychiatric disorders including autism [26-28]. 
A recent transcriptomic analysis of post-mortem brains 
from autistic patients, for example, showed changes in gene 
networks involved in immune and inflammatory responses 
in the frontal and temporal cortices [29]. 



Table 2 Top-ranked biological functions enriched in 
probes differentially expressed (P<0.01) from IPA 



Biological function P value 

Maturation of dendritic cells 8.99 x 1 0"^ 

Influx of phagocytes 1 .69 x 1 0"* 

Binding of interferon-stimulated response element 6.43 x 10"'' 

Damage of oligodendrocytes 7.00 x lO"'' 

Lack of mesenteric lymph node 7.00 x lO"'' 

Lack of peripheral lymph node 7.00 x 1 0"* 

Stimulation of hyaluronic acid 7.00 x 1 0"* 

Morphology of tooth 7.56x10"'' 

Activation of myeloid cells 7.78 x 1 0"'' 

Clustering of lymph node cells 1 .16 x 10"^ 

Osteoclastogenesis of bone cell lines 1.16 x 10"^ 

Transmigration of Thi cells 1.16x10"^ 

Influx of neutrophils 1.21x10"^ 

Proliferation of dendritic cells 1 .27 x 1 0"^ 

Destruction of joint 1.47x10"^ 

Arrest in cell cycle progression of endothelial cells 1.73 x 10"^ 

Lack of cervical lymph node 1.73x10"^ 

Organogenesis of lymphoid organ 1 .73 x 1 0"^ 

Response of carcinoma cell lines 1 .73 x 1 0"^ 

Extravasation of myeloid cells 1 .92 x 1 0"^ 

Morphogenesis of neurites 1 .97 x 1 0"^ 

Induction of neuroglia 2.40x10"^ 

Quantity of IL-5 In blood 2.40 x 10"^ 

Quantity of multinucleated cells 2.40 x 1 0"^ 

Response of lung cancer cell lines 2.40 x 1 0"^ 

Activation of monocytes 2.45 x 1 0"^ 

Acne 2.48x10"^ 

Formation of osteoclasts 2.60 x 10"^ 

Mobilization of blood cells 2.63 x 10"^ 

Activation of phagocytes 2.69 x 1 0"^ 

Binding of lymphoma cell lines 2.76 x 1 0"^ 

Cell viability of blood cells 2.82 x 10"^ 

Quantity of tooth 3.17x10"^ 

Activation of macrophages 3.23 x 10"^ 

Proliferation of B lymphocytes 3.35 x 10"^ 

Response of connective tissue cells 3.42 x 1 0"^ 

Cell movement of dendritic cells 3.88 x 1 0"^ 



In addition to identifying group-level variation be- 
tween the offspring of old and young fathers, we also ex- 
amined family-specific gene expression differences, as 
many of the genetic (or epigenetic) alterations believed 
to underlie the paternal age effect in autism are thought 
likely to be sporadic de novo events occurring in single 
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Figure 2 Top-ranked gene network. The network was derived from Ingenuity Pathway Analysis of differentially expressed genes in the 
offspring of older fathers in a mouse model of advanced paternal age. This network is involved in immune cell trafficking and cell-to-cell 
signalling. Different shapes relate to different molecule types. The key from Ingenuity can be found at [30]. 



litters. We compared the average transcript level for all off- 
spring within each family to the average across all offspring 
of young fathers. We found a 1.8-fold enrichment in the 
number of family-specific significant gene expression differ- 
ences in the offspring of older fathers (average 426 differen- 
tially expressed genes) compared to the offspring of young 
fathers (average 234 differentially expressed genes) {P < 
0.001). Significant probes for each of the six old-father fam- 
ilies are listed in Additional file 3. Of note, a number of 
probes were differentially expressed in the offspring of 
more than one advanced-age father including probes repre- 
senting Mystl, Tnfsfll and Fryl (see Additional file 4). 

Conclusions 

We present evidence for transcriptomic differences in 
the medial prefrontal cortex of offspring of old fathers 
compared to the offspring of young fathers, including 



for genes previously implicated in autism, a neuropsychi- 
atric disease epidemiologically associated with advanced 
paternal age, and an enrichment of loci involved in the 
inflammatory response. Previous studies of gene expres- 
sion changes associated with age in the mouse brain 
have shown enrichment for genes associated with the in- 
flammatory response [8], although this is the first study 
to examine differences associated with advanced paternal 
age. Future work will examine whether these expression 
differences result from de novo genetic (or epigenetic) al- 
terations occurring in the sperm of older fathers. 

Availability of supporting data 

The data set supporting the results is available in the 
Gene Expression Omnibus repository (Currently await- 
ing upload) or downloaded from our webpage http:// 
epigenomicslab.com/Paternal Age data.rar. 
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Additional files 



Additional file 1: Overview of the samples used in this study. 

Additional file 2: Significant gene expression differences for 
individual offspring split by sire for 'young' and 'old' fathers were 
seen for (A) AA408296, (B) MuclS and (C) Lta. 

Additional file 3: Top-ranked differentially expressed transcripts 
(P < 0.001 ) within each of the six families with an old father. Shown 
for each transcript is the corresponding rank in the overall old vs young 
father group comparison. 

Additional file 4: Comparison of the top-ranked differentially 
expressed transcripts (P< 0.001) across each of the six families with 
an old father. Shown for each transcript is the corresponding rank in 
the overall old vs young father group comparison. Empty cells indicate a 
non-significant difference. 



Abbreviations 

CNV: copy number variation; DAVID: Database for Annotation, Visualization 
and Integrated Discovery; FDR: false discovery rate; IL: interleukin; IPA: Ingenuity 
Pathway Analysis; RIN: RNA integrity number; SNP: single nucleotide polymorphism. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

All authors contributed to the design of the study. RS undertook the lab 
work. RS and JM drafted the article. RS collected and analysed data, wrote 
the manuscript, conceived and designed the study, made a critical revision 
and gave final approval of the manuscript CF collected data, conceived and 
designed the study, made a critical revision and gave final approval of the 
manuscript RK collected data, made a critical revision and gave final 
approval of the manuscript. LS analysed the data, made a critical revision 
and gave final approval of the manuscript JB and AR conceived and 
designed the study, provided financial support, made a critical revision and 
gave final approval of the manuscript JM wrote the manuscript, conceived 
and designed the study, made a critical revision and gave final approval of 
the manuscript. All authors read and approved the final manuscript. 

Acknowledgements 

This study was supported by the Beatrice and Samuel A Seaver Foundation, 
by a British Medical Association Margaret Temple Award, and the National 
Institute of Health Research Biomedical Research Centre for Mental Health at 
the South London and Maudsley National Health Service Foundation Trust 
and the Institute of Psychiatry, King's College London, Pilot Award to Drs 
Jonathan Mill and Abraham (Avi) Reichenberg. 

Author details 

'institute of Psychiatry, King's College London, De Crespigny Park, Denmark 
Hill, London SF5 8AF, UK. ^Mount Sinai School of Medicine, Madison Avenue, 
New York, NY 10029, USA. ^University of Fxeter Medical School, St Luke's 
Campus, Magdalen Rd, Fxeter FXl 2LU, UK 

Received: 10 October 2013 Accepted: 10 March 2014 
Published: 23 March 2014 

References 

1. Lundstrom S, Haworth CM, Carlstrom F, Gillberg C, Mill J, Rastam M, 
Hultman CM, Ronald A, Anckarsater H, Plomin R, Lichtenstein P, 
Reichenberg A: Trajectories leading to autism spectrum disorders are 
affected by paternal age: findings from two nationally representative 
twin studies. J Child Psychol Psychiatry 201 0, 51 (7):850-856. 

2. Kong A, Frigge ML, Masson G, Besenbacher S, Sulem P, Magnusson G, 
Gudjonsson SA, Sigurdsson A, Jonasdottir A, Wong WS, Sigurdsson G, 
Walters GB, Steinberg S, Helgason H, Thorleifsson G, Gudbjartsson DF, 
Helgason A, Magnusson OT, Thorsteinsdottir U, Stefansson K: Rate of de 
novo mutations and the importance of father's age to disease risk. 
Nature 2012, 488(741 2):471 -475. 

3. Adkins RM, Thomas F, Tylavsky FA, Krushkal J: Parental ages and levels of DNA 
methylation in the newborn are correlated. BMC Med Genet 201 1, 1 2:47. 



4 Flatscher-Bader T Foldi CJ, Chong S, Whitelaw F, Moser RJ, Burne TH, Fyles 
DW, McGrath JJ: Increased de novo copy number variants in the offspring 
of older males. Trans/ Psychiatr 201 1 , 1 :e34. 

5. Smith RG, Reichenberg A, Kember RL Buxbaum JD, Schalkwyk L Fernandes C, 
Mill J: Advanced paternal age is associated with altered DNA 
methylation at brain-expressed imprinted loci in inbred mice: implications 
for neuropsychiatric disease. Mol Psychiatry 201 3, 1 8(6):635-636. 

6. Colangelo V Schurr J, Ball MJ, Pelaez RP, Bazan NG, Lukiw WJ: Gene 
expression profiling of 12633 genes in Alzheimer hippocampal CA1: 
transcription and neurotrophic factor down-regulation and up- 
regulation of apoptotic and pro-inflammatory signaling. J Neurosci Res 
2002, 70(3)462-473. 

7. Valdes AM, Glass D, Spector TD: Omics technologies and the study of 
human ageing. Nat Rev Genet 2013, 14(9):601-607. 

8. Lee CK, Weindruch R, Prolla TA: Gene-expression profile of the ageing 
brain in mice. Nat Genet 2000, 25(3):294-297. 

9. Paul C, Nagano M, Robaire B: Aging results in differential regulation of 
DNA repair pathways in pachytene spermatocytes in the brown Norway 
rat. Biol Reprod 201 1, 85(6):1269-1278. 

10. Kokkinaki M, Lee TL He Z, Jiang J, Golestaneh N, Hofmann MC, Chan WY 
Dym M: Age affects gene expression in mouse spermatogonial stem/ 
progenitor cells. Reproduction 2010, 139(6):1 01 1-1020. 

1 1 . Alter MD, Kharkar R, Ramsey KF, Craig DW, Melmed RD, Grebe TA, Bay RC, 
Ober-Reynolds S, Kinwan J, Jones JJ, Turner JB, Hen R, Stephan DA: Autism 
and increased paternal age related changes in global levels of gene 
expression regulation. PLoS One 201 1, 6(2):el671 5. 

12. Paxinos G: The Mouse Brain In Stereotaxic Coordinates London: Academic 
press; 2004. 

13. Spijker S: Dissection of rodent brain regions. In Neuroproteomics. New 
York, NY USA: Humana Press Inc Springer; 2011:13-26. 

1 4. Johnson WF, Li C, Rabinovic A: Adjusting batch effects in microarray expression 
data using empirical Bayes methods. Biostatlstics 2007, 8(1):118-127. 

15. Du P, Kibbe WA, Lin SM: Lumi: a pipeline for processing lllumina 
microarray. Bioinformatlcs 2008, 24(1 3):1 547-1 548. 

16. Buyske S, Williams TA, Mars AF, Stenroos FS, Ming SX, Wang R, Sreenath M, 
Factura MF, Reddy C, Lambert GH, Johnson WG: Analysis of case-parent 
trios at a locus with a deletion allele: association of GSTM1 with autism. 
BMC Genet 2006, 7:8. 

1 7. Neale BM, Kou Y, Liu L Ma'ayan A, Samocha KF, Sabo A, Lin CF, Stevens C, 
Wang LS, Makarov V, Polak P, Yoon S, Maguire J, Crawford FL Campbell NG, 
Geller FT, Valladares 0, Schafer C, Liu H, Zhao T, Cai G, Lihm J, Dannenfelser 
R, Jabado 0, Peralta Z, Nagaswamy U, Muzny D, Reid JG, Newsham I, Wu 
YMJ, et al Patterns and rates of exonic de novo mutations in autism 
spectrum disorders. Nature 2012, 485(7397):242-245. 

18. Hu VW, Frank BC, Heine S, Lee NH, Quackenbush J: Gene expression 
profiling of lymphoblastoid cell lines from monozygotic twins discordant 
in severity of autism reveals differential regulation of neurologically 
relevant genes. BMC Genomics 2006, 7:1 18. 

19. Stone JL Merriman B, Cantor RM, Geschwind DH, Nelson SF: High density 
SNP association study of a major autism linkage region on chromosome 
^7.Hum Mol Genet 2007, 16(6):704-71 5. 

20. Al-Ayadhi LY: Relationship between Sonic hedgehog protein, brain-derived 
neurotrophic factor and oxidative stress in autism spectrum disorders. 
Neurochem Res 2012, 37(2):394-400. 

21 . O'Roak BJ, Vives L, Girirajan S, Karakoc F, Krumm N, Coe BP, Levy R, Ko A, 
Lee C, Smith JD, Turner FH, Stanaway IB, Vernot B, Malig M, Baker C, Reilly B, 
Akey JM, Borenstein F, Rieder MJ, Nickerson DA, Bernier R, Shendure J, 
Fichler FF: Sporadic autism exomes reveal a highly interconnected 
protein network of de novo mutations. Nature 2012, 485(7397):246-250. 

22. Hu VW, Nguyen A, Kim KS, Steinberg MF, Sarachana J, Scully MA Soldin SJ, 
Luu T Lee NH: Gene expression profiling of lymphoblasts from autistic 
and nonaffected sib pairs: altered pathways in neuronal development 
and steroid biosynthesis. PLoS One 2009, 4(6):e5775. 

23. Ingenuity Pathway Analysis, [http://www.ingenuity.com/products/ipa] 

24. Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large 
gene lists using DAVID bioinformatlcs resources. Nat Protoc 2009 4(l):44-57 

25. Huang DW, Sherman BT, Lempicki RA: Bioinformatlcs enrichment tools: 
paths toward the comprehensive functional analysis of large gene lists. 
Nucleic Acids Res 2009, 37(1):1-13. 

26. Onore C, Careaga M, Ashwood P: The role of immune dysfunction in the 
pathophysiology of autism. Brain Behav Immun 2012, 26(3):383-392. 



Smith et al. Molecular Autism 2014, 5:24 
http://www.molecularautism.eom/content/5/1/24 



Page 7 of 7 



27. Noriega DB, Savelkoul HF: Immune dysregulation in autism spectrum 
disorder. Eur J Pediatr 20]i, 173(l):33-43. 

28. Micliel M, Schmidt MJ, Mirnics K: Immune system gene dysregulation in 
autism and schizoptirenia. Dev Neurobiol 2012, 72(10):1 277-1 287. 

29. Voineagu I, Wang X, Jolnnston P, Lowe JK, Tian Y, Horvath S, Mill J, Cantor 
RM, Blencowe BJ, Geschwind DH: Transcriptomic analysis of autistic brain 
reveals convergent molecular pathology. Nature 201 1, 474(7351)380-384. 

30. Ingenuity Pathway Analysis legend. [http://w/ww.biolreprod.org/content/ 
suppl/201 l/02/23/biolreprod,l 10.090019.DC1/90019SupLegend.pdfl 

doi:1 0.1 1 86/2040-2392-5-24 

Cite this article as: Smith et at: Transcriptomic changes in the frontal 
cortex associated with paternal age. Molecular Autism 2014 5:24. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at \ rant,,\ 

www.biomedcentral.com/submit Biomea eencrai 



